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The 5N formula for the primordial curvature perturbation C, is extended to include vector as well 
as scalar fields. Formulas for the tree-level contributions to the spectrum and bispectrum of C, are 
given, exhibiting statistical anisotropy. The one-loop contribution to the spectrum of C, is also worked 
out. We then consider the generation of vector field perturbations from the vacuum, including the 
longitudinal component that will be present if there is no gauge invariance. Finally, the 5N formula 
is applied to the vector curvaton and vector inflation models with the tensor perturbation also 
evaluated in the latter case. 
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I. INTRODUCTION 

Starting at an 'initial' temperature of a few MeV, the observable Universe is now understood in considerable detail. 
At the initial epoch the expanding Universe is an almost isotropic and homogeneous gas. The perturbations away 
from perfect isotropy and homogeneity at the initial epoch are the subject of intense study at present, because they 
determine the subsequent evolution of all cosmological perturbations According to observation, the dominant and 
perhaps the only initial perturbation is the curvature perturbation so-called because it is related to the perturbation 
in the intrinsic curvature of space-time slices with uniform energy density. 

To understand the nature and origin of C, one uses comoving coordinates x, that move with expansion of the 
unperturbed Universe. Also, one considers the Fourier components with comoving wave-vector k. Physical positions 
are a(i)x and physical wave- vectors are k/a(i), where a is the scale factor of the Universe. The Hubble parameter is 
H = a/a, with a dot denoting derivative with respect to the cosmic time t. 

It is convenient to smooth all relevant quantities on a comoving scale, somewhat below the shortest scale of 
cosmological interest. This will not affect the Fourier components on cosmological scales, and will greatly simplify 
the analysis. Consider a given cosmological scale, characterised by wavenumber k/a. On the assumption that gravity 
slows down the expansion of the cosmic fluid, aH/k = d/k increases as we go back in time. At the present epoch 
scales of cosmological interest correspond to 10~^ < aH/k < 1, but at the 'initial' temperature T ~ MeV they all 
correspond to aH/k 1. Such scales are said to be outside the horizon. 

To explain the origin of the perturbations, it is supposed that going further back in time we reach an era of inflation 
when by deflnition gravity is repulsive. At the begining of inflation the smoothing scale is supposed to be inside the 
horizon. With mild assumptions, it can be shown that inflation drives all perturbations to zero at the classical level. 
But as each scale k leaves the horizon, the quantum fluctuations of those scalar fleld perturbations with mass m < H 
are converted 0, @| to classical perturbations. 

According to the usual assumption, one or more of these scalar fleld perturbations is responsible for the curvature 
perturbation (for a recent account with references see Ref. [3|). In that case, the statistical properties of C (specifled 
by its correlators) are homogeneous and isotropic (invariant under displacements and rotations). It has been pointed 
out recently that vector fleld perturbations could contribute to ( 0, 0,0, Such contributions will typically make 

C statistically anisotropic, but still statistically homogeneous. 

It was shown in an earlier pape r 11(1 h ow, including only scalar flelds, one may calculate the correlators of C through 
what is called the SN formalism [ill, ■ The 6N formalism has recently been applied to the vector fleld case in 
a particular setup In this paper, we work out a completely general SN formalism including vector flelds and then 
apply it to a different setup used for the vector curvaton iS;l2i| and vector inflation [l3| scenarios. 
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The plan of the paper is the following. In Section HIl we give some useful formulas and survey the observational 
status regarding statistical anisotropy. Section lUll is devoted to a brief description of the 6N formaHsm, this time 
including vector fields. In Section |IV] we calculate the spectrum of C at tree and one-loop level, and the bispectrum 
of C at tree level. In Section |V] we recall the generation of a scalar field perturbation from the vacuum. In Section IVTl 
we see how a gauge field perturbation can be generated. In Section fVIII we see how a vector field perturbation can 
be generated, using a modified-gravity action without gauge invariance and including the longitudinal component. 
In Sections IVIIII and [IX] we see how a vector field perturbation can contribute to C, through respectively the vector 
curvaton and vector infiation mechanisms. We conclude in Section [Xl 



II. OBSERVATIONAL CONSTRAINTS ON THE CURVATURE PERTURBATION 

Direct information on the curvature perturbation comes mostly from measurements of the anisotropy of the CMB 
and the inhomogeneity of the galaxy distribution. These cover a limited range of scales, corresponding to roughly 
Aln/c ~ 10 where k is the comoving wavenumber. Indirect information is available at much longer and shorter scales. 
In this section we summarise the information. 



A. Formulas 



We are interested in the correlators of the curvature perturbation, in particular the two-point correlator. For any 
cosmological perturbation /?(x), at some fixed time, we define Fourier components with normalisation 



/?(k) = J /3(x)e-*'^ ''d'^a; . (1) 

Assuming that the two-point correlator (/3(x)/3(x')) is invariant under translations (statistically homogeneous), the 
two-point correlator of the Fourier components takes the form 

27r2 

(/3(k)/3(k')) = {2nf6{k + k')^PMk) , (2) 

which defines the spectrum If the two-point correlator is also invariant under rotations (statistical isotropy) 

the spectrum 7^/3 (k) depends only on the magnitude k. In that case we shall sometimes invoke a quantity P^ik) = 
{27ry k^)Pp{k). 

By virtue of the reality condition /3(— k) = (3*{k), an equivalent definition of the spectrum is 

(/3(k)/3*(k')) - {27Tf6{k - k')^Vp{k) . (3) 

Setting k = k' the left hand side is (|/3(k)p). It follows that the the spectrum is positive and nonzero. 

Even if Pfl(k) is anisotropic, the reality condition requires Vpik) = Pf3{—k). The anisotropy will therefore be of 
the form (Tsj (see also Refs. [lB,[i3|) 

T'Mk) = 7'i°(fc) [l + g^(d.k)2 + ...] , (4) 

where 'Pp°{k) is the average over all directions, d is some unit vector and k is a unit vector along k. 

If there is no correlation between the Fourier components except for the reality condition, the perturbation is said 
to be Gaussian. Then the two-point correlator is given by Eq. ^ and the three-point correlator vanishes while the 
four-point correlator is 

(/3ki/3k./?k3/?k4> = (/3ki/3k.)(/3k3/3k4) + (/3ki/3k3>(/?k./3k4> + (/3ki/3k4> (/3k./3k3> • (5) 

The five-point correlator vanishes and the six-point correlator is given by the analogue of Eq. ([5]), and so on. All 
correlators are known once the spectrum is specified. We conclude that a Gaussian perturbation is statistically 
homogeneous even though it need not be statistically isotropic. 



*^ The averages are over some ensemble of universes, of which our observable Universe is supposed to be a typical realization. 
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Non-gaussianity is signalled by a non-vanishing 3-point correlator, an additional ('connected') contribution to the 
4-point correlator and so on. Statistical homogeneity requires that each correlator of Fourier components vanishes 
unless the sum of the wave- vectors vanishes (generaHsing the delta function of Eq. Q), and statistical isotropy requires 
that it is invariant under rotations. In particular, statistical homogeneity requires a 3-point correlator of the form 

(/3(k)/3(k')/3(k")) = {2Trf5{k + k' + k")B(k, k', k") , (6) 

and statistical isotropy requires that B depends only on the magnitudes of the vectors. Assuming statistical isotropy 
one also defines a reduced bispectrum Bp by 

Bpik, k', k") = Bp{k, k', k") [P[}{k)Pp{k') + cyclic permutations] . (7) 



B. Spectrum and non-gaussianity 

Observational results concerning the spectrum Vq are generally obtained with the assumption of statistical isotropy, 
but they would not be greatly affected by the inclusion of anisotropy at the 10% level. 

Direct observation, coming from the anisotropy of the CMB and the inhomogeneity of the galaxy distribution, gives 
information on what are called cosmological scales These correspond to a range A In fc ~ 10 or so downwards 

from the scale k~^ ~ ^cT^ that corresponds to the size of the observable Universe^^. It is found that is almost 
scale independent with the value Vq^^'^ ~ 5 x 10~^. There is mild scale dependence corresponding to 

n - 1 = ^^^^ -0.040 ± 0.014 . (8) 
din K 

On much bigger or smaller scales the constraint is far weaker. Assuming a constant n on such scales, they are 

50 

- 5 < n - 1 < 0.4— , iVcorr = ln(fccorr/fcmax) ■ (9) 

*corr 

The lower bound, referring to very large scales k <^ Hq, comes [l| from the absence of an enhancement of the CMB 
quadrupole (Grishchuk-Zeldovich effect). 

The upper bound is more interesting. In this expression. Neon is the number of e-folds of inflation, between horizon 
exit for the smallest cosmological scale fc";^^ a-nd horizon exit for the smallest scale k'^^^ on which the curvature 
perturbation exists (correlation length). It corresponds [l^ to the following values for the spectrum at those scales: 

n'^'(fcmax) < 5 X 10-5 , 7'c'/'(fccorr) < 10"^ . (10) 

The first number is the observed value on cosmological scales. The second number corresponds to an order of 
magnitude upper bound on the spectrum that under certain assumptions is required to avoid an overabundance of 
primordial black holes ^d\. Further discussion about the upper bound on is given in Ref. [l3|- 

If C is generated during inflation, or soon afterwards, kcon will be the scale leaving the horizon at the end of 
inflation. Then A'corr — — 10, where A^ is the number of e-folds of inflation after the largest cosmological scale H^^ 
leaves the horizon. For a high inflation scale and a fairly standard cosmology afterwards, A^ ~ 60 making A^corr — 50. 
If instead ( is formed long after inflation, through say the curvaton model, A'corr can be much lower for the same A^, 
and A^ itself will be reduced if the inflation scale is low. 

If the spectral tilt varies, the upper bound refers to average of the tilt with respect to In fc, in the interval fcmax < fccorr- 
The possibility of large tilt on small scales has been investigated in Ref. A strongly increasing tilt on small scales 
could come from a single mechanism for generating n, such as the running mass inflation model. Alternatively, a large 
and practically constant n on small scales could be generated if the curvature perturbation has two components: 

n(^) =^flat(fc)+^stccp(fc). (11) 

The flrst component might be nearly flat and dominate on cosmological scales, while the second might have large tilt 
and dominate in the interval fcmax < k < fccorr- In that case, the upper bound in Eq. |[9)) appHes to the spectral tilt 

of T^stccp- 



As usual a subscript indicates the present epoch, and we set ao = 1. 
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Coming to non-gaussianity, one generally focusses on the bispectrum, working with the quantity /nl = (5/6)S(;. If 
/nl is generated from one or more gaussian field perturbations with scale-independent spectra it is practically scale 
independent. With that assumption, the most recent analysis [2l| finds /nl = 38 ± 21 at Icr but —4 < /nl < 80 
is allowed at 95% confidence level. For fully correlated non-gaussianity, /nl'Pc"^^^ of order the fractional non- 
gaussianity of C which means that the non-gaussian fraction is less than 10~^ or so, and in any case the observational 
bound on /nl corresponds to a small non-gaussian fraction [j^ . 

Allowing scale dependence of the bispectrum, the observational bounds are very weak on scales outside the cosmo- 
logical range, so that for example C. could be the square of a gaussian quantity. 

C. Statistical anisotropy and statistical inhomogeneity 

Taking all the uncertainties into account, observation is consistent with statistical anisotropy and statistical in- 
homogeneity but allows either of these things at around the 10% level. In this section we briefly review what is 
known. 

Assuming statistical homogeneity of the curvature perturbation, a recent study (see also Refs. (23 . [25| l of 
the cosmic microwave background radiation (CMB) temperature perturbation flnds weak evidence for statistical 
anisotropy. They keep only the leading term of Eq. (0]) : 

7'c(k)=7'f°(fc)(l+g(d.kf) , (12) 

and flnd g ~ 0.15 ± 0.04 with d in a specifled direction. The authors point out though that systematic uncertainties 
could make g compatible with zero. We will therefore just assume \g\ < 0.3 In other words, we assume that the 
spectrum of the curvature perturbation is isotropic to within thirty percent or so. There is at present no bound on 
statistical anisotropy of the 3-point or higher correlators. 

In some different studies, the mean-square CMB perturbation ir i oppos ite hemispheres has been measured, to see 
if there is any difference between hemispheres. A recent work [13, Hail^l flnds a difference of order ten percent, for 
a certain choice of the hemispheres, with statistical signiflcance at the 99% level. Given the difficulty of handling 
systematic uncertainties it would be premature to regard the evidence for this hemispherical anisotropy as completely 
overwhelming. 

Let us see what hemispherical anisotropy would imply for the curvature perturbation. Focussing on a small patch 
of sky, the statistical anisotropy of the curvature perturbation implies that the mean-square temperature perturbation 
within a given small patch wih in general depend on the direction of that patch. This is because the mean square 
within such a patch depends (in the sudden decoupling approximation) upon the mean square of the curvature 
perturbation in a small planar region of space perpendicular to the line of sight located at last scattering*^. But 
the mean-square temperature will be the same in patches at opposite directions in the sky, because they explore the 
curvature perturbation ((k) in the same k-plane and the spectrum V(^{'k) is invariant under the change k — k. It 
follows that statistical anisotropy of the curvature perturbation cannot by itself generate a hemispherical anisotropy. 

In the above discussion of the CMB temperature perturbation, we ignored cosmic variance, by identifying the 
measured mean-square temperature perturbation within a given patch with the ensemble average of that quantity. 
That will certainly be permissible if the multipoles of the CMB, including the lowest ones, are almost uncorrelated 
corresponding to an almost gaussian curvature perturbation. 

With the caveat concerning cosmic variance, we conclude that hemispherical anisotropy of the CMB temperature 
requires statistical inhomogeneity of the curvature perturbation. Then (C(k)C(k')) is not proportional to (5(k-|-k'). But 
in a small region of the observable Universe it might still be reasonable to invoke approximate statistical homogeneity, 
by defining a position-dependent spectrum V{k,x) (taken for simplicity to be rotationaiiy invariant). This way of 
generating the hemispherical anisotropy has been considered in Refs. [30l . [3l| . but is outside the framework of the 
present paper. 

Before ending this section we note that, in addition to the primordial curvature perturbation, there might be a 
primordial tensor perturbation with spectrum Vh jH. The fraction r = Vh/VQ is constrained by observation to be 
<0.1 0. 



A related work [2^ shows that the lowest detectable value for |g| from the expected performance of WMAP is \g\ ~ 0.1. The same 
analysis gives the lowest detectable value from the expected performance of PLANCK: \g\ ~ 0.02. 
*^ The sudden decoupling is not essential here. It can be replaced by the exact line of sight formalism, leading to the same conclusion. 
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III. THE 5N FORMALISM 

The 6N formalism for scalar field perturbations was given at the linear level in Refs. (ill , [l^ . At the non-linear 
level which generates non-gaussianity it was described in Refs. [l3, Here we extend the formalism to include 
vector fields. 

With generic coordinates the line element of the perturbed universe is 

ds^ = Of^^dx^'dx" . (13) 

The coordinate system of the perturbed universe defines a slicing (constant time coordinates) and a threading (constant 
space coordinates) of spacetime. 

To define the cosmological perturbations, one chooses a coordinate system in the perturbed universe, and then 
compares that universe with an unperturbed one. The unperturbed universe is taken to be homogeneous, and is 
usuahy taken to be isotropic as well. In this Section though, we develop the SN formahsm without assuming isotropy. 

The SN formalism does not invoke a theory of gravity, but it does invoke an energy- momentum tensor T^i^ . From a 
mathematical viewpoint, any definition will do provided that it satisfies the continuity equation S/ ^iTj^ — with Vm 
the covariant derivative. Following for instance Refs. H^l, we define T^i/ in terms of the spacetime curvature: 

R^u - ^9,.uR = -SnGTf,, . (14) 

This is the Einstein field equation if, in a locally inertial frame, is the energy-momentum tensor of Special 
Relativity. In the context of field theory, this means that the action should be of the form 



ImlR + L 



(15) 



where mp = (SttG)"^/^ is the reduced Planck mass, and C, evaluated in a locally inertial frame, is the lagrangian 
density of fiat spacetime field theory. Then T^^, is the 'improved energy-momentum tensor' which is given in terms of 
the fields by a standard expression. The bosonic part £bos of L gives a contribution 

dC 

^^'^°' = 2^^-5p>''Cbos. (16) 

Of course, we can always write the action in the form given by Eq. (fT5| with some L. When that is done, the 
contribution of the bosonic part £bos will still be given by Eq. iflGj) . We shall invoke this expression in several cases 
where Einstein gravity holds, and will invoke it in Section [IX] for a case where Einstein gravity does not hold, dropping 
the label 'bos'. 



A. The curvature perturbation and the tensor perturbation 

To define the curvature perturbation, we smooth the metric tensor and the energy-momentum tensor on a comoving 
scale significantly shorter than the scales of interest, and we consider the super-horizon regime aH ^ k. On the 
reasonable assumption that the smoothing scale is the biggest relevant scale, spatial gradients of the smoothed metric 
and energy-momentum tensors will be negligible. As a result, the evolution of these quantities at each comoving 
location will be that of some homogeneous 'separate universe'. In contrast with earlier works on the separation 
universe assumption, we will in this section allow the possibility that the separate universes are anisotropic even 
though homogoneous. 

We consider the slicing of spacetime with uniform energy density, and the threading which moves with the expansion 
(comoving threading). By virtue of the separate universe assumption, the threading will be orthogonal to the slicing. 
The spatial metric can then be written as 

5,,(x,T)^a2(x,T)(/e2'^(-'-)) , (17) 

where / is the unit matrix, and the matrix h is traceless, which means that /e^'* has unit determinant. The time 
dependence of the locally defined scale factor a(x, t) defines the rate at which an infinitesimal comoving volume V 
expands: V/V = 3a/ a. 
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We split In a and hij into an unperturbed part plus a perturbation: 

In a(x,T) = In a(T) + C(x,r) , (18) 
hij{x,T) EE hij{T) + Shij{x,T) . (19) 

The unperturbed parts can be defined as spatial averages within the observable Universe, but any definition will do as 
long as it makes the perturbations small within the observable Universe. If they are small enough, C and Shij can be 
treated as first-order perturbations. That is expected to be the case, with the proviso that a second-order treatment 
of C will be necessary to handle its non-gaussianity if that is present at a level corresponding to /nl ^ 1 (with the 
gaussian and non-gaussian components correlated) [3^ . 



1. The curvature perturbation 

In this paper we are mainly concerned with the curvature perturbation C Because le^ has unit determinant, 
the energy continuity equation d{Vp) = —PdV implies that C, is independent of position, during any era when the 
pressure P is a unique function of the energy density p [l3| (hence uniform on sHces of uniform p) . Absorbing ^ into 
the unperturbed scale factor, C(x) is then time independent. 

From the success of Big Bang Nucleosy thesis, we know that Einstein gravity is a good approximation when the 
shortest cosmological scale approaches horizon entry at T ~ 1 MeV. Also, the cosmic fiuid is then radiation dominated 
to high accuracy implying P — p/3 and a constant value of We denote this value simply by C(x), and it is the one 
constrained by observation as described in Section Ull 



2. The tensor perturbation 

The perturbation Shij may also be of interest. We discuss it at this point in general terms, and in Section HXl we 
provide an explicit calculation within the vector infiation model. 

Consider first the unperturbed quantity hij{T). In this paper we are taking the unperturbed expansion to be 
practically isotropic expansion with Cartesian coordinates. As a result, we can take the unperturbed quantity to vanish 
so that a(T) is the unperturbed scale factor. More generally, if the unperturbed quantity is any time-independent 
matrix, we can make a Hnear coordinate transformation which diagonaHses le'^ and can then choose the normalization 
of the scale factor so that hij again vanishes. A time-dependent unperturbed quantiy hij[T) would correspond to an 
unperturbed Universe with anisotropic expansion. 

If one or more vector fields exist during infiation, one might think that the unperturbed expansion may e asily be 
anisotropic. Assuming Einstein gravity though, that is not the case because according to a theorem of Wald [3J, |33| 
enough infiation driven by a constant scalar field potential will isotropise the expansion*''. This statement becomes 
only an approximation for realistic slow roll infiation where the potential is varying, and it doesn't apply to 'vector 
infiation' models where infiation is driven by a constant vector field potential [14, 36, 33, [s^l- For vector infiation 
though, one can ensure approximate isotropy of the expansion by invoking a large number of independent fields [l^ . 
as we shall discuss in Section ITXl 

After infiation, an era of anisotropic stress (from vector fields or any other source) might cause significant anisotropy 
of the expansion, but that does not happen in the usual scenarios of the early Universe. Assuming Einstein gravity, 
the anisotropy will anyway decay when the anisotropic stress switches off. 

As we are deahng with a smoothed metric well after horizon exit, the status of the perturbed quantity /iij(x, t) 
at a given location is the same as that of the unperturbed quantity. The anisotropy of the local expansion will be 
neghgible if Wald's theorem holds or if there is vector infiation with a sufficiently large number of independent vector 
fields. Then the perturbation 6hij{x,T) is almost time independent. Also, we expect Shij to remain time-independent 
after infiation until the approach of horizon entry, since the local anisotropic stress is expected to remain negligible. 

Now we consider first order cosmological perturbation theory, taking the unperturbed hij to vanish. At first order, 
the equations satisfied by the cosmological perturbations comprise three uncoupled modes, termed scalar, vector and 
tensor. The first order perturbation Sgij is equal to Sij( + Shij with ( belonging to the scalar mode. Setting spatial 



It is so-called because one usually has in mind the case that Shij is negligible; of course it too corresponds to a perturbation in the 
spatial curvature. 

He calls this constant potential a cosmological constant. 
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gradients equal to zero in accordance with the separate universe assumption, Shij satisfies the transversality condition 
diShij = 0, which means that it belongs to the tensor mode. It can be written in terms of polarization tensors as 
etjh+ + e~jh-, and assuming statistical parity invariance each of the amplitudes has the same spectrum 7^tcn/4. The 
fraction r = Vtcn/Vc; is < 10^^ and future measurements will reduce this bound by a factor of 10 to 100, or detect 
r [si. 

Let us discuss the origin of the tensor perturbation Shij, within first order cosmological perturbation theory assuming 
Einstein gravity. The standard calculation assumes negligible anisotropy in the infiationary expansion and negligible 
anisotropic stress (which will certainly be the case if scalar fields dominate). Under these assumptions, each of 
hj^, x /V2mp has the action of a free scalar field. The classical equation of motion is 

6h,j + SHShij + {k^la^)5hij = 0, (20) 

which makes 5hij constant after horizon exit. The spectrum of the perturbation generated from the vacuum fiuctuation 
is 

which is too small to observe in small-field infiation models. 

If one or more vector fields is relevant during infiation, the unperturbed expansion will be anisotropic at some 
level. Then, even if the vector field is unperturbed, 5hij will time dependent during infiation, and can be generated 
even if the vacuum fiuctuation is negligible; such a contribution would be correlated with the curvature perturbation 
[13, SH, in contrast with the one generated from the vacuum fiuctuation. A perturbation of the vector field will give 
an additional effect. As we argued earher these effects are expected to be neligible if Einstein gravity holds, assuming 
that either Wald's theorem applies or there is vector infiation with sufficiently many independent fields. 

Within this first order treatment, the tensor perturbation is gaussian. Since the tensor perturbation has yet to 
be detected there is little motivation to consider its non-gaussianity. At the time of writing, the only calculation 
of non-gaussianity has been done by Maldacena [i^l assuming single field slow roll infiation with Einstein gravity. 
Using second order perturbation theory he chooses a gauge where 5hij is transverse as well as traceless. He calculates 
the three-point correlators involving Fourier components of C, and/or at the epoch soon after horizon exit, and 
finds them to be suppressed by slow roll factors. If C, receives contributions only from the infiaton perturbation, it 
is constant after horizon exit and then the three point correlator of C, corresponds to /nl ^ 10~^ which is almost 
certainly too small ever to detect. There is no reason to think that the correlators involving 5hij will be detectable 
either. Judging by this example, there is no need for the discussion of Shij to go beyond first order cosmological 
perturbation theory. 



B. The SN formula 



Keeping the comoving threading, we can write the analogue of Eq. ifTT)) for a generic sHcing: 

5,,(x,r)^a2(x,r)(/e2''(-'-)) , (22) 

with again /e^'' having unit determinant so that the rate of volume expansion is V/V — 3a(x, t). Starting with an 
initial 'fiat' slicing such that the locally-defined scale factor is homogeneous, and ending with a slicing of uniform 
density, we then have 

C(x,t) =(57V(x,t), (23) 

where the number of e-folds of expansion is defined in terms of the volume expansion by the usual expression N = 
V/3V. The choice of the initial epoch has no effect on i5iV, because the expansion going from one flat slice to another 
is uniform. We will choose the initial epoch to be a few Hubble times after the smoothing scale leaves the horizon 
during inflation. According to the usual assumption, the evolution of the local expansion rate is determined by the 
initial values of one or more of the perturbed scalar flelds 4>i . Then we can write 

<^/(x) = 0/ + J07(x), (24) 
C(x,t) = 57V(0i(x),02(x),...,i) 

= Ni{t)64>i{^) + \Nij{t)5M^)5<t>j{^) + ..., (25) 
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where Nj = dN/d(f>i, etc., and the partial derivatives are evaluated with the fields at their unperturbed values denoted 
simply by 0/. The field perturbations (50/ in Eq. l(25|) are defined on the 'fiat' slicing such that a(x, t) is uniform. 

The unperturbed field values are defined as the spatial averages, over a comoving box within which the perturbations 
are defined. The box size aL should satisfy LHq ^ 1 so that the observable Universe should fit comfortably inside it 
[43 |. If there have been exponentially many e-folds of infiation before the observable Universe leaves the horizon, one 
could choose In(LHo) to be exponentially large, but that would not be a good idea because it introduces unknowable 
new physics and places the calculation out of control [43|. One therefore chooses a 'minimal box', such that ln{LHo) 
is significantly bigger than 1 without being exponentially large. 

The spatial averages of the scalar fields, that determine Nj, etc., and hence C cannot in general be calculated. 
Instead they are parameters, that have to be specified along with the relevant parameters of the action before the 
correlators of C can be calculated. The only exception is when C is determined by the perturbation of the infiaton in 
single-field infiation. Then, the unperturbed field value when cosmological scales leave the horizon can be calculated, 
knowing the number of e-folds to the end of infiation which is determined by the evolution of the scale factor after 
infiation. Although the unperturbed field values cannot be calculated, their mean square for a random location of the 
minimal box (ie. of the observable Universe) can sometimes be calculated using the stochastic formalism 

In this paper we suppose that one or more perturbed vector fields also affect the evolution of the local expansion 
rate. Keeping for simplicity one scalar field and one vector field we have 

C(x, t) = <5^(0(x), Ai^),t) = N^6(b + N\6A, + ^N^^i6(b)^ + dA, + ^N'JJA, SA, + ... , (26) 

where 

Ar - . _ dN _ d^N .J _ d^N . _ d^N ^ ^ 

with i denoting the spatial indices running from 1 to 3. As with the scalar fields, the unperturbed vector field values 
are defined as averages within the chosen box. 

In these formulas there is no need to define the basis (triad) for the components Ai. Also, we need not assume that 
Ai comes from a 4-vector field, still less from a gauge field. 

The discussion so far allows the unperturbed expansion to be anisotropic. In the following sections though, we will 
take it to be isotropic. Also, we take the unperturbed spatial geometry to be fiat. Then the unperturbed line element 
is 

ds^ = a2(r) {-dr^ + Sijdx'dx') , (28) 

where r is conformal time and a is the scale factor. Depending on the context, we may instead use cosmic time t 
corresponding to dt = adr. We shall take Ai to be the physical field, defined with respect to the orthonormal basis 
induced by the Cartesian space coordinates = a{t)x''. We shall also have occasion to consider the field Bi — aAi 
that is defined with respect to the orthogonal (but not orthonormal) basis induced by the comoving coordinates x\ 
The corresponding upper-index quantities are A^ — Ai and ~ a^^Bi. 



C. The growth of C 

As noted earlier, ( is constant during any era when pressure P is a unique function of energy density p. In the 
simplest scenario, the field whose perturbation generates C is the infiaton field (p in a single-field model. Then the 
local value of is supposed to determine the subsequent evolution of both pressure and energy density, making ( 
constant from the beginning. 

Alternatives to the simplest scenario generate all or part of C. at successively later eras. Such generation is possible 
during any era, unless there is sufficiently complete matter domination {P — 0) or radiation domination {p = P/S) 
Possibilities in chronological order include generation during (i) multi-field infiation (ii) at the end of infiation 
(iii) during preheating, (iv) at reheating, and (v) at a second reheating through the curvaton mechanism (i^. [ZtI. l48l. 

A vector field cannot replace the scalar field in the simplest scenario, because unperturbed infiation with a single 
unperturbed vector field will be very anisotropic and so will be the resulting curvature perturbation. Even with 
isotropic infiation, we are about to see that a single vector field perturbation cannot be responsible for the entire 
curvature perturbation (at least in the scenarios that we discuss) because its contribution is highly anisotropic. It 
could instead be responsible for part of the curvature perturbation, through any of the mechanisms listed above. Of 
these, the end of infiation mechanism has already been explored In this paper we explore another one, namely the 
vector curvaton mechanism [5]. We will also explore the vector infiation scenario |i.l4], according to which infiation is 
driven by a large number of randomly oriented vector fields which can give sufficiently isotropic infiation and (as we 
shall see) an extremely isotropic C. 
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IV. FORMULAS FOR THE SPECTRUM AND BISPECTRUM OF THE CURVATURE 

PERTURBATION 

A. Spectrum of the vector field perturbation 

In Section |V] we describe the standard scenario for generating the scalar field perturbations from the vacuum. 
Within this scenario, these perturbations are Gaussian with no correlation between different perturbations. Their 
stochastic properties are defined by the spectrum Vstj} of each field. Either of the equivalent definitions ^ and ^ 
can be used to define the spectrum, with (3 — 6(f). 

To deal with a vector field perturbation SAi we write 

6 A, (k, r) EE ^ (k)MA (k, r) , (29) 

A 

where with the z axis along k the polarization vectors are defined by 

= (l,z,0)/\/2, e« = (l,-^,0)/\/2, e^^s = (0, 0, 1) . (30) 

These expressions define the polarization vectors only up to a rotation about the k direction but that is enough for 
the present purpose. We will let the change k — k reverse z and x but not y. Then eA(~k) — — e^(k) and there is 
a reality condition v4^(k, r) — — Aa(— k, t). 

If the vector field corresponds to a gauge field, we choose the gauge so that ^long — leaving only Al and Ar. 
Otherwise we have to keep all three ^a- 

In Sections I VII and rVIII we describe two scenarios for generating the vector field perturbations SA\. Within both of 
them, these perturbations are statistically isotropic and Gaussian, with no correlation between different A or between 
the perturbations of different fields (scalar or vector) . As a result we need only to consider the spectra Vx = VsAx ■ 
They can be defined by the analogue of either Eq. ([2]) or Eq. 

(MA(k)MKk')) = (27r)3^(k-k')|^7'A(fc), (31) 

(MA(k)MA(k')) = -{2Tif5{\<i + \<^)'^Vx{k). (32) 

The spectra are nonzero and positive, with the minus sign in the second expression coming from eA(— k) = — e^(k). 

We will normally have Vl = Vr, since a difference between these quantities would indicate parity violation of the 
evolution of Ai . It is therefore useful to define 

V± = ^iVR±VL) , (33) 

so that only V+ will normally be present*^. 

In the models that we discuss, the scale dependence of the spectra Vxik) comes from the evolution of the perturbation 
5 Ax after horizon exit during infiation. In this regime, the spatial gradient k/a is negligible compared with the Hubble 
parameter, and we expect that it will be neghgible compared with any other relevant parameter*^. In that case, the 
evolution of (5Aa(x, r) at each position will be the same as for the unperturbed field ^i(r). By rotational invariance 
the evolution of the latter is independent of i. Therefore, we expect that the evolution of the three perturbations 
6Ax will become the same after horizon exit, giving them the same spectral index. In that case Hong, defined as 
Hong = 'Piong/'P+, will be just a number, independent of k. 

The correlators of the 6Ai (k) are 

(M,(k)M,(k')) = (27r)3<5(k + k0|f f7^™'^(k)P+(fc) +^7^°<^<^(k)P_(A:) +T;™^(k)7'io„g(fc)l ' (34) 



where 



i;7™(k) = <5y - fc,^„ T°j'^ik) = e,,kkk, Tjp^ik) = kkj. (35) 



Calculations that generate V- as well are described in Refs. |5ol.[5H. 
This is verified for the specific scenarios that we consider. 
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B. Spectrum of 

1. Tree-level spectrum 

Since C is gaussian to high accuracy, it seems reasonable to expect that C will be dominated by one or more of the 
linear terms in Eq. (|26| . Keeping only them (corresponding to what is called the tree-level contribution) we find*^° 

T'f °<=(k) = NlVsm + N\N\ [7-™Xk)7'+(fc) + T^"S(k)7'io„g(fc)] (36) 
= NlVs^ik) + NlV+{k) + (N^ • ]^fV+{k) (no„g - 1) . (37) 

The above corresponds to Eq. (fT2l) with d = N^i, N^i being the Cartesian vector with components N\, and 

Vr{k) = N^Vs^k) + NlV+ik) , (38) 



NlV+{k) 
N^Vs^ik) + NlV+(k) ' 



(nong - 1) n^::lL n.. > (39) 



where Na = \/ ^a^a magnitude of N^. The spectrum is scale- invariant if the spectra of the field perturbations 

are scale invariant. 

If the vector field perturbation dominates C we have simply g = riong — 1. If the vector field is a gauge field Hong = 0, 
and if its action is Eq. (|78l) below r-jong = 2. In both cases, the the observational bound \g\ < 0.3 is violated which 
means that the vector field contribution cannot dominate. If there is no other vector field contribution, the dominant 
contribution to C must then come from one or more scalar field perturbations. 

To avoid the need for scalar perturbations, one can suppose that a large number N of vector fields perturbations 
contribute to C, with random orientation of the unperturbed fields. With a sufficient number of fields, there is then 
no preferred direction and the curvature perturbation is isotropic. 



2. One-loop contribution 
Using Eq. |[5]), the contribution from the quadratic terms (one- loop contribution) is 

+ \nIaKaPs4>(\^ + p|) [7--(p)7'+(p) + T^"S(p)7'iong(p)^ 
+\NtAKA{nr{^ + p)T;r"(p)^+(|k + v\)V+{p) + 
+T°,dd(k + p)T°i^\^)VA\^ + v\)V-{v) + 

+ p)T;r^(p)Plo„g(|k + pDT'longlp) + 

+2rr"(k + p)7;'r (P)n(|k + p|)7'io„g(p)}} . (40) 

If the spectra are scale-independent, the integral is proportional to In(fcL) [53] where L is the box size. If we 
allow ln(A:L) to be exponentially large the one-loop contribution can dominate the tree-level contribution even with 
C almost gaussian, but the whole calculation is then out of control [i^. With a 'minimal' box size such that In(fcL) 
is not exponentially large, and keeping only a single scalar field contribution, it has been shown that the ratio 
(P^i-ioop/P^trcc)i/2 -g ^j^g fractional non-gaussianity /nl'T^c^^^ ^he curvature perturbation which from 

observation is < 10"'^. However, the loop contribution to C from a given field could dominate the tree level from that 
field, if both contributions are small compared with the total. This could in particular be the case for the vector field 
contribution. 



The terminology tree-level and one-loop corresponds to a Feynman graph formalism [s^ l that could easily be extended to include vector 
fields. 
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C. Bispectrum of 



Working to leading order in the quadratic terms of the SN formula, we arrive at the tree-level contribution to the 
bispectrum. Evaluating it using Eq. ^ we find 



Bf^"(k,k',k") = NlN^4Ps^{k)Ps^{k')+ eye. perm.] + 



5 perm. | 



X "7-r"(k')p+(fc') + d(k')p-(fc') + r;r'(k')^iong(fc'; 

where Ps4,{k) and P\{k) are defined as 



eye. perm 



}• 



1^ 



(41) 



(42) 



Reversal of the three wave- vectors corresponds to the parity transformation, and from the reality condition C(^k) = 
C* (k) it changes each correlator into its complex conjugate. For the spectrum this is not of interest because the reality 
condition also makes the spectrum real. For the bispectrum with statistical isotropy it is also not of interest, because 
the reality condition plus statistical isotropy make the bispectrum real*^^. In our case, the bispectrum is statistically 
anisotropic, and is guaranteed to be real only if the parity-violating spectrum V- vanishes. 

Existing analysis of the bispectrum assumes statistical isotropy*^^, and it seems important that the analysis should 
be extended to allow for anisotropy and possible parity violation. The relation between non-gaussianity and the 
anisotropy of the spectrum is explored in Ref . [H^ . 

The second order contribution of the quadratic terms in the SN formula gives the one-loop contribution to the 
bispectrum. It could be significant or even dominant. It has been calculated for the scalar case in Ref. and has 
been investigated for the case of multifield infiation in for instance Refs. [H^, [H3|- The one- loop contribution from a 
vector perturbation will be given in a separate pubHcation [H3|- 



V. SCALAR FIELD PERTURBATION FROM THE VACUUM FLUCTUATION 

During infiation, both scalar and vector field perturbations can be generated from the vacuum fiuctuation. We 
begin by describing carefully the scalar field calculation, emphasising some points that will be important when we 
come to the vector field. 



A. General considerations 



We shall focus on the simplest setup. Only the few e-folds either side of horizon exit are considered. Unperturbed 
infiation is supposed to be isotropic, and almost exponential so that the Hubble parameter can be taken to be 
constant. It is assumed that the field perturbations can be treated as free fields, so that they satisfy uncoupled linear 
field equations. Also, the scalar fields are taken to five in unperturbed spacetime, which means that the back-reaction 
of the fields on the metric is ignored. By virtue of these features, the scalar field perturbations are gaussian and 
statistically independent, and the object of the calculation is to calculate their spectra. 

For the calculation itself we do not need to invoke a theory of gravity or a model of infiation. But these things are 
needed if one wishes to check that the back-reaction is neghgible and the field is practically free. Assuming Einstein 
gravity and slow-roll infiation, the check has been done as follows. First, the modification of the linear evolution 
equation to include back-reaction has been calculated, both for single-field [5§| and multi-field infiation. It is 
found to be small, provided that the relevant fields are slowly varying on the Hubble timescale, as will be the case if 
their potential is fiat enough for the slow-roll approximation to apply*^^. Second, the treatment of the perturbation 



The triangle of vectors obtained by reversing the vectors can be brought into coincidence with the original triangle by a rotation. 
See for instance Ref. [H. 

From the form of the back-reaction, one expects this to be the case even for non-Einstein gravity 
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has been carried out to second order [i^ . IgiI . Ig^ and third order [g^ . [63 | (including the back- reaction) . From this the 
3-point and connected 4-point correlators of ( were calculated. They were found to be negligible in accordance with 
the linearity assumption. 

These calculations invoke only scalar fields, which is consistent with the assumption of isotropic unperturbed 
infiation. In the present paper we are going to suppose that one or more vector fields exist during infiation. As we 
noticed in Section IIII A 21 the unperturbed (spatially homogeneous) part of a vector field will at some level cause 
anisotropic unperturbed expansion. This will break the rotational invariance of the evolution equations for the scalar 
[15, 40, 41, 65, 66J, causing their spectra to be anisotropic. At the moment it is not understood how to calculate 
the spectra of scalar field perturbations in such a case, because the linear evolution equations have singlular solutions 
[gsI . [6^ . The generation of vector field perturbations will also be affected by anisotropic unperturbed expansion, 
though that has yet to be investigated. As we saw in Section IIII A 2^ the level of anisotropy in the expansion is 
expected to be small and in this paper we simply ignore it. 



B. Quantum field theory 

There is no need to assume Einstein gravity during infiation. We need only the effective action for the scalar 
field, vaHd while relevant scales are leaving the horizon. We describe the standard scenario, in which is canonically 
normalized. Although the calculation works for a more general potential, it will be enough here to consider the 
quadratic case: 



The action is then 



S" = i J dTd^Xy/^[Ccf,{T,x) 



(43) 

(44) 
(45) 



This is supposed to hold to good accuracy while scales of interest leave the horizon, with practically constant 
during that era. The dots indicate contributions, which generate infiation if that is not already done by (j) 

As the field is supposed to live in unperturbed spacetime described by the fine element in Eq. (|28)) we can write 



-d^<pd''(l>-^a\T)m^ 



(46) 



with the index fi now raised by ry^'^ instead of g^^" . 
The unperturbed field equation is 



(47) 



where an overdot denotes d/dt. We take infiation to be practically exponential so that a oc exp{Ht). We assume that 
(p is a light field, defined as one with 



(48) 

—m'^<j)/3H, which is expected to hold more or 
less independently of any initial condition. Then the fractional change in (f) over one Hubble time is much less than 
1. If the inequality is only marginally satisfied it will be of order 1. 
For the first order perturbation we work with ip = aS(f>. It satisfies 



If this inequality is well satisfied there will be a slow roll solution 



ip"{k, t) + (fc2 + a^m^) <p(k, r) = 



m 



2H^ 



(49) 



If it is done by <f> there is slow-roll inflation with 4> the inflaton, but we are not assuming slow-roll inflation and still less that </> is the 
inflaton within that paradigm. 
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where a prime denotes d/dr. To arrive at the quantum theory we need the action for 54>, obtained from Eq. (|45l) . 
After dropping a total derivative it is 



(50) 



drd^k 



^'^{\i,T)-{e + a'm^)^{\i,T) . (51) 



For each k this is the action of an oscillator with time-dependent frequency. 

We adopt the Heisenberg picture whereby the state vector is time independent. Promoting Lp to an operator we 
write 

^[a(k)^(fc,r)e'''- + at(k)^*(fc,r)e-*-] . (52) 

The mode functions 'p{k, t) satisfy the same evolution equations as the classical perturbations (/^(k, t). The former are 
independent of the direction of k because the evolution equations do not pick out a preferred direction, and neither 
does the initial condition that we come to shortly. 

The consistent quantization of this system requires the commutation relation 

[a(k),at(k')] = (2^)3j(k-k'), (53) 

and the Wronskian 

(p*{k,T)drip{k,T) - ip>{k,T)dr(p*{k,T) = -i. (54) 
Well before horizon exit, ip is a linear combination of exp(±ifcr). We make the usual choice 

(55) 



which will be justified shortly. We postulate a unique vacuum state, annihilated by the a(k), and take the Hilbert space 
to be Fock space, whose basis is built by acting on the vacuum by products of the creation operators d^(k). The basis 
vectors are eigenvectors of the occupation number operator = L^'^aj^Sk, which gives the number of particles with 
momentum k. The particle interpretation can be justified using Eq. lfT6|) for the energy momentum tensor. It shows 
that the vacuum state, with zero occupation number, has momentum density and pressure pvac = —Pvac = A'^/IGtt^ 
where A is the ultra-violet cutoff. This is set equal to zero by absorbing it into the scalar field potential. Then the 
energy- momentum tensor of a generic basis state is that of a gas of particles with the relevant occupation numbers. 
If the occupation numbers depend only on the direction of k, the momentum density and anisotropic stress vanish, 
leaving pressure and energy density P — p/3. 

The final step is to assume that the time-independent state vector is close to the vacuum state. In other words, we 
assume that the the occupation number Uk of the quantum states (averaged over a cell of k space) is much less than 1. 
With Einstein gravity, that assumption is mandatory if there have been AA^ ^ ln(Mp/7J,) e-folds of infiation before 
cosmological scales leave the horizon, because the positive pressure P ^ nk{k/a)'^ from particles with momentum of 
order k/a would otherwise overwhelm the negative pressure P = —3MpH^ that is required for infiation P, [g^]. The 
condition AA^ >• ln(Mp/i/*) is quite mild, and will almost certainly be satisfied for the shortest cosmological scale if 
infiation takes place at the usual high scale iJ* ~ 10~^Mp. 

Instead of using the negative frequency mode function in Eq. l(55|) , one might consider using a linear combination 
of positive and negative frequencies. This corresponds to using annihilation operators dk, related to the original ones 
by a BogoHubov transformation: 

dk = afcdk -I- PkoL I (56) 

with |afcp = 1-1- |/3fcp. A Fock space vector, labelled by the eigenvalues of rik — al^ak/L^, does not have well-defined 
nk and does not have well-defined energy-momentum tensor either. In a state where rik has expectation value (?ik), 
the expectation value of nk is 



("k) = (^k) + |/3fe|' (l + 2(7^k)) 



(57) 
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The expectation value of the energy-momentum tensor in this state is that of a gas with occupation number (rik) . As 
in the previous paragraph, it is reasonable to require this occupation number is much less than 1, in order to ensure 
that the positive pressure of the gas will not be significant at the beginning of inflation. Looking at Eq. (|57|) . we see 
that this requires |/3fe| ^ 1. In words, the initial mode function cannot be much different from the negative frequency 
mode function in Eq. fSS]) "^^^. 

This argument for the choice of the negative frequency solution relies on the fact that it minimizes the energy density 
and the pressure, of the gas of particles that will be present if any other choice is made. The standard argument 
[g^ I invokes only the energy density, which by itself would not be dangerous. Indeed, one is already discounting the 
vacuum energy density Pvac, which is permissible because it comes with Pvac = —Pvac- 



C. Spectrum of the perturbation 

To calculate the spectrum of (p, we identify the ensemble average in Eq. ([2]) as a vacuum expectation value, with (p 
replaced by (p. Then 



fc3 



r^ik,T) = \cp{k,T)\' 



(58) 



Apart from the reality condition, there is no correlation between different Fourier components, because there is no 
correlation between their vacuum fluctuations and no coupling between their evolution equations. In other words, the 
perturbation ip is Gaussian in the linear approximation that we are using. 

The mode function is the solution of Eq. (|49|) with the initial condition Eq. l(55|) . For m = it is 



(pik,T) 



i (kr — i) 
2k kr 



(59) 



Well after horizon exit this gives 



27r J 



(60) 



Keeping to, we can write Eq. l(49|) as 



di- 



(^(k,r) =0, 



(61) 



With#16 



9 / 



H 



(62) 



This is the Bessel equation with independent solutions Ju{kT) and J-^{kT). The solution satisfying the initial 
condition is 



TT e 2 



(63) 



Well after horizon exit this gives 



' ^- 23/2r(|)V2fc^ ' 



(64) 



The field becomes classical in the sense that [pik, r), i9r<^(k, t)] tends to zero provided that v is real. This condition 
corresponds to m? < jH^ . We reject the regime m? <C —H^ because the spectrum is too steep to be of interest. 



Of course, it also requires that the state is close to the vacuum state, corresponding to (n^) <t. 1. 
As indicated we choose the positive sign. 
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In any case, the calculation almost certainly becomes invalid in this regime for two reasons. First, the unperturbed 
field (f) will roll rapidly away from the origin, making it unHkely that the neglected terms of the potential in Eq. l(43|) 
remain neghgible over the several Hubble times that it takes for relevant scales to leave the horizon. Second, the 
back-reaction of the perturbation on the metric will probably not be negligible. These are the considerations that 
require the light field condition in Eq. l(48|) . Applied to negative m^, this condition is equivalent to v^l. 



Well after horizon exit Eq. (|64|) gives 



{l-cos2-Kv) \2ttJ \2aH 



2tt ) \aH 
2m 



(65) 



2 



"-scalar - 1 ^ ^ '^'^ ~ SIP ' ^^^^ 

The final equality is valid for |m^| ^ H^. 

Instead of taking 4> to be massless, one might think that a more accurate early-time approximation would be 
obtained by keeping the mass m. That is not the case though, because the classicality condition in Eq. (|48|) means 
that the effect of m is no bigger than the effect of the expansion rate H. Because m is so small, we cannot regard it 
as a mass term in a fiat spacetime quantum field theory. 



D. Spectral tilt 

The scale dependence given by Eqs. i(65l) and l(66|) can be understood in the following way. Soon after horizon exit, 
when the classical perturbation first emerges, its spectrum is roughly independent of m, and hence ~ {H/2'k)'^. After 
that, the perturbation evolves according to Eq. l(49|) with k^ = 0, which gives the second factor of Eq. (|65l) . 

Even more simply, we can understand the scale dependence just from the unperturbed equation l|47p . It is the same 
as Eq. l(49l) with k — and has two independent solutions. One is proportional to a2(--i) and the other to a ' " '. 
The second solution decays relative to the first by a factor a"'*'^, and one expects that it will become negligible soon 
after horizon exit*^^. Using the first solution we again arrive at Eqs. (|65|) and l(66|) . 

The curvature perturbation is given in terms of the scalar field perturbations by Eq. (|25l) . Let us take the initial 
epoch in that equation to be after all cosmological scales have left the horizon, but not too long after. Then the 
estimate Eq. l(65l) should apply to each scalar field perturbation. Supposing that a single scalar field perturbation 
dominates and using the tree- level expression for , the spectral index n of C will obviously be equal to the spectral 
index Ugcaiar of the scalar field. The observed spectral tilt value n — 1 ~ —0.04 suggests that the light field condition 
in Eq. ijiS]) is very well satisfied by the relevant field. 

Since the tree-level expression for 7^^ treats the field perturbations Hnearly, one can instead calculate the spectral 
index of C using the 'horizon-crossing trick', whereby the initial epoch is instead taken to be a fixed number of Hubble 
times after horizon exit for the scale k. This technique allows one to easily include a slow variation of H, defined 
by en = —H/H^. It reduces n by an amount 6€h if (p is the infiaton and by 2eH otherwise. The horizon crossing 
technique also allows one to write down a formula for n if several scalar fields contribute, in terms of the first and 
second derivatives of the potential at horizon exit P, 



VI. GAUGE FIELD PERTURBATION FROM A TIME-DEPENDENT GAUGE COUPLING 



In this section and the next, we see how a vector field perturbation may be generated. In this section we work with 
the following effective action during almost-exponential infiation: 



S 



(67) 



'^^'^ Unless V is close to 1 corresponding to "scalar — 4. 
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where F^^, = dfj^B^, — d^B^ is the field strength with a gauge field. It can be written 

S-- 



drd^a 



(68) 



where now the indices are raised with 77''^ instead of g^'^ . 

If / is time-independent it can be set equal to 1 because any constant value can be absorbed into B^. Otherwise, / 
represents a time-dependent gauge coupling. To respect invariance under time displacement, / should be a function 
of one or more fields with no explicit time dependence. 

As with the scalar field, there is no need to assume Einstein gravity during infiation. The other terms in the action 
are supposed to give infiation with practically constant H, and to generate /(r) without having any other effect on 
the evolution of the gauge field during infiation. For that to be the case, any scalar field coupled to must have 
zero value (no spontaneous symmetry breaking) with negligible quantum fiuctuation around that value. 

Starting with Ref. j6§|, this action has been widely considered for the generation of a primordial magnetic field, and 
it has recently been considered Q for the generation of a vector field perturbation that can generate a contribution 
to (. In the latter context, an extension to include a mass term is studied in Refs. [5, 70]. 

By a choice of gauge we set Bq and djB^ equal to zero. We assume almost exponential infiation and work with the 
perturbation 



A = fSB^ = aSA, 



(69) 



We are absorbing / into the definition of the physical field Ai even though it is supposed to be varying while 
cosmological scales are leaving the horizon. At some stage / will become time-independent making Ai indeed the 
physical gauge field. 

The perturbation has only transverse components, which satisfy the field equation 



r 
f 



A(k,r) =0, 



(70) 



with X = L or R. The prime denotes d/dr. 

The quantization is just like the scalar case [tH. Each A\ has the scalar field action in Eq. (|5T1) . with {aifif' 
replaced by —f"/f- We write 



A(x,r)= / (0jE[^'Wa,(k)A(fc,r)e 



+ ef (k)al(k)^*(fc,T)e- 



with the sum going only over \ — L,R. The commutator is 

aA(k),4,(k')l = (27r)35(k - k')<5AA' 



(71) 



(72) 



and the Wronskan of A\{k,T) is —i. Well before horizon exit /"// is supposed to be negligible and one adopts the 
initial condition 



-ikr 



2k 



(73) 



as well as the Fock space, and one assumes that the state is close to the vacuum state. These assumptions can be 
justified in the same way as for the scalar field case. 

Following Refs. [1, oj we adopt the parameterisation / oc a". Then Eq. l(70|) has the same form as Eq. l(6T1) for the 
scalar field perturbation: 



9^ 



^A(fc,r) = 0, 



with v = 



Well after horizon exit, it leads to a classical perturbation with the spectrum 

fc3 1 



Vx{k,T) 



27r2 a2 



\Ax{k,T)\\ 



(74) 



(75) 
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Using the solution of Eq. (|74l) with the initial condition in Eq. (|73|) we have 

1 



(76) 
(77) 



The spectrum is scale invariant if a = —2 or a = 1 

Since v is always real, the vacuum fluctuation always gives a classical perturbation after horizon exit. We reject 
V ^ \ (equivalent to a ^ 1) because the predicted spectrum is too steep to be of interest. 

As was pointed out in Ref. [t^I, a classical perturbation is obtained even with the standard gauge coupling 
corresponding to a = 0. In that case the evolution of the mode function is not affected by horizon exit and nvcc — 1 = 2. 
This can be traced to the fact that the action is invariant under a conformal transformation of the metric, which 
means that we can go to the flat spacetime metric. After horizon entry during the post-inflation era, classicality is 
lost and we recover the vacuum state of the late-time quantum field theory, but that is of no concern in the present 
context. Of course it prevents one using the standard action to generate a primordial magnetic field (quite apart from 
the fact that the spectral index would anyway be too big for the field to be useful). 

Taking H to be constant, the contribution of the vector field contribution to Q has spectral index index rivoc- As in 
Eq. (flTI) the vector contribution could dominate on small scales, and even the conformal invariant tilt rivoc — 1 = 2 
might be allowed by the bound in Eq. ^ though that would need a rather low value Nik^aa-x) ~ 10. 



VII. VECTOR FIELD PERTURBATION WITH COUPLING TO R 

A. The action 

As an alternative to the previous case, we now consider the following effective action during infiation: 



S = / dTd 



J3„ 



(78) 



The third term of this action violates gauge invariance. As a result, one cannot use gauge invariance to motivate 
the particular form of the kinetic term, and one cannot use any other internal symmetry either. The most general 
quadratic kinetic term consistent with Lorentz invariance is [lH| 

/:kin = -piV^B'^V^B, - /?2 (V^B^)' - PsV^B'^V.B^ , (79) 

with V being the covariant derivative. Gauge invariance requires /3i = — /33, which is the only restriction provided by 
symmetry considerations. The action in Eq. I|78p invokes that condition, without the justification of gauge invariance. 

The motivation for the action in Eq. l(78|) comes, not from symmetry considerations but because it has two 
remarkable properties. One property concerns the perturbation 5B^ that is generated from the vacuum fiuctuation. 
As we will show in this section, the spectrum of the perturbation is scale-invariant if m = 0, for both the transverse 
and longitudinal perturbations. This calculation of the spectrum invokes no theory of gravity. The other remarkable 
property concerns the theory of gravity and will be described in Section |IX] (generalizing the action to include an 
arbitrary number of vector fields) . These special properties perhaps suggest that the action in Eq. l(78|) can emerge 
in a natural way, in the context of field theory or perhaps string theory. 

Much of the literature, starting with Ref. [73|, goes further and identifies the field i?^ in Eq. iffS]) with the 
electromagnetic field. That requires its couplings to other fields (including the known Standard Model fields) to be of 
the standard gauge-invariant form even though there is no gauge invariance*^^. It seems to us to be a step too far, 
when one can as well generate a primordial magnetic field using the gauge invariant action of the previous section. 

We require the other terms of the action to generate infiation, without affecting the evolution of B^ during infiation. 
For that to be the case, any terms coupling B^ to scalar fields should have a negligible effect. There is no reason to 



In Ref. 01 this is given incorrectly as a = —1. Note that the value a = 2, advocated in Ref. |7ll | in the context of a primordial magnetic 
field, makes the energy density rather than the field perturbation scale invariant. 

The form of the coupling of the photon to spin half fields is completely determined by renormalizability, but not the form of its coupling 
to the and Higgs fields. 
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suppose that such coupHng occurs through the gauge-invariant terms of the form —T>^4>{'D^(f>)*. But if for instance a 
(global or gauge) U{1) symmetry acts on the phase of (j) but not on i?^ one might have a term of the form — |0pi?^_B^ 
and then we are requiring that the the C/(l) is unbroken with neghgible quantum fluctuation, just as in the gauge- 
invariant case except that the U{1) now has nothing to do with B^. 



B. Generating the field perturbation 

As the action in Eq. (|78| contains no time derivative for the time component Bq, this component is related to the 
space components Bi by a constraint equation*^". We take the spacetime metric to be unperturbed. 

The unperturbed field has zero time component, and the space components of the physical field Ai — Bi/a satisfy 

A, + 3HA, + m^A,=0. (80) 

This is the same as for a scalar field with mass-squared m^. 

As in the previous section, we work with the perturbation of the physical field, Ai = aSAi = SBi. We expand its 
operator in the form given by Eq. ifflj) . including now the longitudinal mode since there is no gauge invariance. 

Consider first the transverse modes, X = L,R. They satisfy the equation 0,01 

[d^ + a^m^ + k^]Ax=0, (81) 

where*^^ 

rh'^ =m'^ + ^R = - 2H'^ . (82) 

This is the same as for a scalar field with mass-squared . The action for each of A\ is also the same [6^ . We adopt 
the initial condition, the Fock space, and the vacuum state assumption, with the same justification as in the scalar 
field case. Then 



nvec-1 = 3-2^^— , (84) 

A classical perturbation is generated if v is real corresponding to rr? < 9H^/4. As with the scalar case, we reject the 
case ^ —H^. The spectrum is too steep to be of interest, and anyway the evolution of Ai would be so rapid that 
additional terms in Eq. (fTS)) (required to stabiHze Ai) could hardly remain neghgible over the several Hubble times 
that it takes for cosmological scales to leave the horizon. We therefore require 

-H'<m'< . (85) 

As advertised, the tilt vanishes if to = 0. 

Now we discuss the quantization of the longitudinal perturbation. Its mode function satisfies 0| 



Aong = . (86) 

For m = corresponding to = —2H^, the independent solutions (given here for the first time) are 



Aong(fc^) « [-kr + ^ ± 2*) e^^'^^ . (87) 
We see that the solutions are regular even at the point where the round bracket in Eq. l(86|) vanishes. 



#20 pqj. generic choice of the kinetic term, Bo becomes an independent field. Its perturbation is considered in Ref. [zliIzHl- 
We used the relation -R = —12H'^, valid during exponential infiation. 
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We can show that the solution of Eq. (|86l) is non-singular even for ^ 0. This can be done by using the Frobenius 
method for differential equations with regular singular points (see for example Ref. [Z^). First we make a change of 
variables 



y 



a m 



1 



with y varying in the region —l<y< oo. Eq. l(86|) with this transformation translates into the form 



" 2y(y+l) ^ 



1 ~ 9 1 






y 



4(y + lf 



(89) 



with < and the regular singular point at y ^ 0. The general solution of this equation can be found using the 
ansatz 



A 



long 



n=Q 



Dny 



s+n 



(90) 



where Dq ^ 0. In this case the series in Eq. l(90l) is convergent at least in the region — 1 < y < 1 without a singular 
point. We will show that it converges even at this point and that the ansatz in Eq. I|90p gives two independent 
solutions. To show this let us substitute Eq. ((90| into Eq. (|89p giving 



n=0 



Dn 



A{s + n) {s + n - 2)y' 



s+n-2 



+ 8{s + n)[s + n--] y'+'''^ + 



+4 (s + n) ( s + n - I ) y'^" 



s+n+l 



0. 



(91) 



In order for the equality in Eq. ((911) to be valid, coefficients in front of each y with the same power must vanish. 
The coefficient in front of the term with the smallest power, i.e. y"*"^, is 4£>os (s — 2). Because ^ 0, from the 
indicial equation s (s — 2) = we find 



s = , or s — 2. 



(92) 



Because these two solutions differ by an integer, it might be alarming that the general solution of Eq. (|89|) might 
involve the logarithm. However, by closer inspection of Eq. (|9T|) we find that the coefficient D2 of the series with 
s = is arbitrary, thus the power series in Eq. l(90|) with s = and s = 2 give two independent solutions. And 
because the series does not involve negative powers of y, i.e. s > 0, it converges at the singular point y ~> 0. 

The action corresponding to Eq. l(86|l is*^^ 



i / drd^kC, 



C ~ [arhY' 



l-^;ong(k,r)|^ 



|Aong(k, r) 



To set the initial condition well before horizon entry we define A — {a\fn\/k) Along- In the regime a|TO| <C fc. 



C^±{\AY -k'\A\' 



(93) 
(94) 

(95) 



where the sign ± is that of fh? , hence negative for the case of interest ~ —2H^. 

Except for the negative sign this is same as for the scalar field case. To quantize it we assume the same initial 
condition A = exp^—ikr) / V2k , and adopt the vacuum state. The justification for these assumptions is similar to the 



This is given for the case m? = —2H^ in Ref. |66l |. and it can be derived by perturbing the full action. Of course it is unique only up 
to a total derivative. 
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one that holds for the scalar field (and transverse vector field), but not identical because of the negative sign. Because 
of this sign, occupied initial states would have negative energy density and pressure, P = p/3 ~ — ^^(/c/a)^. As the 
pressure is negative it is not dangerous for inflation. Instead, it is the negative energy density that is dangerous. As 
the total energy density is required to be positive, the negative contribution of occupied states has to be less than the 
total at the beginning of inflation. Assuming as before AiV ^ Mp/if* e-folds of inflation before cosmological scales 
leave the horizon, this again requires occupation number much less than 1, justifying both the choice of initial mode 
function and the assumption of the vacuum state. 

The spectrum T'long is given by Eq. ((75l) . For m = 0, corresponding to rh^ = —2H^, we flnd well after horizon exit 



This corresponds to Hong = 2, which according to the discussion at the end of Section flV B II means that the vector 
fleld perturbation cannot generate the dominant contribution to the curvature perturbation. 

It has been suggested [sl, [6^ that the action in Eq. (|94l) does not correspond to a well defined quantum field theory 
for negative ifi^. We have demonstrated that there is a well defined quantum field theory even in this case. Before 
the epoch |amp = k^, a negative rfi^ corresponds to a negative kinetic term in the action. This will cause some 
degree of instability when more terms are included in the action, corresponding to the interaction of Aong with other 
fields and/or gravity. But such interactions are assumed to be negligible whenever one considers the generation of a 
gaussian classical field perturbation from the vacuum fiuctuation, and as we mentioned already has been justified for 
both scalar and vector field perturbations. In this connection, it is important to realise that the the negative sign 
holds only before the epoch |amp = which is around the time of horizon exit. Also, that only a limited number of 
e-folds of infiation take place between the emergence of k/a from the Planck scale and horizon exit, which means that 
there is only a limited amount of time for the presumably small interactions of -4iong to have any effect. After horizon 
exit, the evolution at each location is given by the classical expression in Eq. ((80l) and we have no more need of the 
quantum theory. According to the classical expression Ai is slowly varying. It moves towards zero if is positive. 
If instead is negative moves towards the vev of Ai. That vev will be at the minimum of the potential V{Bf^B'^), 
whose leading term rn^B^B^^ /2 is displayed in the action in Eq. l(78|). 



We have described two mechanisms that can generate a vector field perturbation from the vacuum fiuctuation. 
In this section and the next we describe two mechanisms by which such a perturbation can give a contribution to 
the curvatu re p erturbation. We begin in this section with the vector curvaton mechanism [11 . This is the curvaton 
mechanism (S, 53, H, , using a vector field instead of the usual scalar field. 

The vector curvaton field Ai (x, r) is smoothed on a scale somewhat below the shortest cosmological scale and it 
has a perturbation a~^Ai = 5Ai. After horizon exit during infiation, the spatial gradient of Ai becomes negligible 
and it evolves at each point as an unperturbed field. In the simplest curvaton scenario, which we adopt, the evolution 
is negligible during and after infiation, until some epoch when Ai begins to oscillate. At this epoch, there is supposed 
to be Einstein gravity and the effective action is supposed to be 



This is the action of a massive vector field, living in the expanding Universe which is taken to be unperturbed. When 
H falls below the mass m, the field begins to oscillate with angular frequency m. As the spatial gradient is negligible, 
the oscillation is a standing wave whose initial amplitude varies with position. 

As originally proposed, the vector curvaton scenario generates the perturbation 5Ai with essentially the action in 
Eq. (|78l) . taking to be a constant parameter which during infiation is negligible. For the present purpose there is 
no need to say how the perturbation is generated. 

The energy density of the oscillation is, in terms of the physical field Ai — Bi/a, 




(96) 



VIII. VECTOR CURVATON 




(97) 




(98) 



(99) 
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where Ostart is the scale factor just before the oscillation starts. In the second line, A is the unperturbed value just 
before the oscillation starts and (5A(x) is its perturbation. The oscillation amplitude falls Hke a~^/^, and ispractically 
constant during one oscillation. As a result, the stress is practically zero just as in the scalar field case We take 
the decay to be instantaneous, which from the scalar field case we know will be an adequate approximation. 

The contribution of pA to the total energy density is supposed to be initially negligible, and with it the contribution 
C,A of 5Ai to C- But the oscillation is supposed to take place in a radiation background, so that Pa/ P grows Hke a{t) 
and Ca becomes significant. 

To calculate C,a we will use the following expression f49|: 

a = W-^. (100) 

3 PA 

_ _3pA_ PA , . 

dpA + 4pr P 

where p = pA + Pr- This expression is vaHd to first order in SpA, which is evaluated on a 'fiat' slice where a(x, t) is 
unperturbed. 

We take the curvaton to decay instantly (sudden-decay approximation) and evaluate (a just before the curvaton 
decays, assuming that C is constant thereafer. The final equality in Eq. I|10ip is justified because the sudden decay 
approximation gives an error of similar magnitude, both errors disappearing in the Hmit ^Ia = 1- Evaluating 5pA to 
first order we have 

Ca^^^a^. (102) 



3 " 



The tree-level contribution to the spectrum is 



4 n 



^C.(fc) = 9^^+(fc) [1 + (^long - 1) (A • k)^J , (103) 
where A = A/|yl|. 

The spectum V+ik) is to be evaluated just before the oscillation starts. In Ref. [l3| it is taken to be the same as 
that at the initial epoch during infiation and that in turn is supposed to be generated from the action in Eq. ((78|. 
Then 7^+ is given by Eq. (|83l) with rivcc practically equal to 1. 

Evaluating SpA to second order we have [l3| 



2^ A6A, , 1 6 AS A, ^^^^^ 



This is valid only for ft a <C 1. To handle the case Qa — 1 one could go to second order in 5pA, or much more simply 
evaluate N and hence SN directly*^^. All of this is the same as for a scalar field contribution, where the evaluation 
of N was done in Ref. [l3|. We shall not pursue the case Qa — 1 in the present paper. 

Our Eq. lfT02l) is Eq. (64) of Ref. 0, generalized to allow Ua < 1 and written to exhibit manifest invariance under 
rotations. The spectrum V^^ was not calculated in Ref. but it was implicitely assumed to be rotationally invariant 
so that it could be the dominant contribution. 

In accordance with the discussion at the end of Section IIV B 1\ this realisation of the vector curvaton mechanism 
cannot give the dominant contribution to (. It could do so by invoking^ several vector curvaton fields. We note that 
the case of several scalar curvaton fields has been considered in Ref. 

IX. VECTOR INFLATION 

Recently, it has been proposed (see also Refs. [tI, [tI, [sOI) that infiation can be driven by a large number of 
independent vector fields. They considered only the unperturbed case, and invoked the large number to make the 
unperturbed metric practically isotropic. We consider the perturbation. 



#23 To first order in SpA, one finds by that method N\ = 2QyiAi/3|Ap, in agreement with Eqs. H102II and II104II . 
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The action is Eq. l(78|). extended to include many vector fields: 



S 



E 

b 



4 fJ-" 



-r] S^bC')^ 



6 



(105) 



As it is supposed to apply throughout inflation (starting with the approach of horizon exit for the largest cosmological 
scale k ~ Hq), the additional terms are supposed to be negligible throughout that era, and not just while cosmological 
scales are leaving the horizon. Also, the action is supposed to define the theory of gravity as well as the dynamics of 
the vector fields. 

Consider first the unperturbed fields Bf^lr). Because each of them has a direction, the expansion is not generally 
isotropic but the anisotropy can be negligible if there is a large number of randomly oriented fields which is 
assumed. Given a large number of fields, the randomness assumption is well justified because, as stated in Section [TTH 
the unperturbed field values are defined as spatial averages within a chosen box, whose location is random. By the 
same token, it does not seem reasonable to replace the randomness assumption by the assumption that there are three 
fields whose unperturbed values are orthonormal, though that would also give unperturbed spacetime (stI. IssI^^'^. 

Varying the action with respect to an unperturbed field, one finds that Eq. ijSOl) is satisfied. Varying the action 
instead with respect to the spacetime metric gives the right hand side of the Einstein field equation, which we take 
as the definition of the energy momentum tensor. For a generic spacetime, the term coupling R to the vector fields 
would make the form of this energy momentum tensor dependent on the metric; in other words it would modify 
Einstein gravity. Remarkably though, the modification is neghgible when spacetime is practically unperturbed 
As a result we have the usual expressions, depending only on the vector field: 



b,i 



A 



(106) 
(107) 



The Friedmann equation therefore takes the usual form, irripH^ ~ p. 

From Eqs. ((80l) . I|106p . and (|107p we see that each component of the unperturbed field is equivalent to a scalar field. 
In the regime > there is infiation, with 



6 ni p ^-^ 



(108) 



It follows that the number of e-folds to the end of infiation is given by the same expression as in the scalar field case 
0,111: 



N ■ 



Ami ^ 



(6) 1 2 



(109) 



Now we consider the curvature perturbation generated by vector infiation. It turns out to be practically the same as 
if the field components are replaced by scalar fields and that case has already been worked out using the 5N formaHsm 
[sH. The derivatives of N for use in the 5N formula are given by Eq. I|109p : 



A 



(b) 



- 2ml 



1 



2m% 



SijSab 



(110) 



The transverse spectrum of the field perturbations are given by Eqs. l(76|) and (|77l) (the same as for a scalar field) 
and the longitudinal spectra are Piong = 2V+. 

The spectrum is given by Eq. ((57|) (without the scalar contribution) , summed over all of the vector fields using 
Ny^(b) — A^^^/2TOp. Since ^ H^, we have V+ ~ (-fffc/27r)^ for each field, where Hk is the Hubble parameter when 



The choice might be justified on anthropic grounds if isotropic expansion was favoured on those grounds but that there is no suggestion 
that such is the case. In particular there is no suggestion that the 30% or so of anisotropy allowed by present data is anthropically 
disfavoured. 
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the scale k leaves the horizon. Since there are a large number of randomly oriented fields we can pretend that they 
all have the same magnitude when evaluating the second term. Since the average of cos^ is 1/2, this gives 



. (- 



2 



VAk) - -N — ^ . (Ill) 

Except for the factor 3/2, the spectrum is the same as was found for the scalar field case [81]. Such a result is 
independent of the number of fields. 

Assuming that TV ~ 55 e-folds of infiation take place after the observable Universe leaves the horizon, the observed 
magnitude of Vc^ is reproduced ii H ~ 10^** GeV at the end of infiation*^^. The non-gaussianity is neghgible and the 
spectral index isn=l — 2/iV. 

In Ref. fl^ (see also Ref. [sOl) the tensor perturbation Shij is also considered, actually for a wide class of vector 
infiation models including the one considered here. The tensor perturbation is supposed to live in unperturbed 
spacetime, as in the standard calculation described in Section IIII A 21 But because the action in Eq. I|105p does 
not correspond to Einstein gravity its linear evolution equation differs from Eq. (|20|) . It is found that Shij can have 
significant time dependence, which makes it difficult to see how the prediction can be compared with observation. 

To avoid these problems, one should go to the Einstein frame by making a conformal transformation of the metric. 
Suppose that we have an action of the form 



j drd^X^ [fR + /:n.atter(gM-: •••)], (112) 

with / being any scalar function of bosonic fields, and £mattcr a function of the bosonic fields that is obtained from a 
fiat spacetime expression through the replacements 77^1/ — > .9^,^ and 9^ where Vm the covariant derivative*^^. 

Now we make a conformal transformation of the metric, g^v — exp(2il)5p^ with exp(2r2) = 2f /mj^. After dropping 
a total derivative this gives |82] 

drd^xV^f + Anattcr ) (113) 



20 , 



-Cmatter = 2 Vp^^ V^^^ + 6^ Anattcr (ff^"^ (Spi. ) , ' ' •) • (114) 

The usual application (s^ . [83| is to slow roll infiation, with f{(f>) a function of just the infiaton field (j>. Then the 
conformal transformation just gives (j) a non-canonical kinetic term. The single field can be redefined to have a 
canonical kinetic term, so that we again have slow roll infiation though with a different potential. 
In our case, 



b 

Since / is slowly varying, there is almost exponential infiation in the Einstein frame just as in the original frame, with 
practically the same Hubble parameter i?*. Now though, £mattcr is a complicated function, making the Einstein frame 
completely unsuitable for the calculation of the vector field perturbations, and hence of the curvature perturbation. 
It is however the one in which one should calculate the tensor perturbation. 

Indeed, the standard first-order cosmological perturbation theory calculation described in Section IlII A 21 will apply 
in the Einstein frame, provided that the number of vector fields is large enough to make the anisotropic stress tensor 
negligible. The linear evolution of the tensor perturbation in the Einstein frame is then given by Eq. (|20l) . and the 
spectrum is given by Eq. I|2ip . The tensor fraction r = T^ton/T'c therefore given by the same formula as in the scalar 
field case, which is |8l| r — 8/N. 

Unfortunately these combined predictions for n and r are disfavoured by observation [l3|. Making the masses 
unequal would make the spectral index even less than one without altering r [8l|, which increases the disagreement 
with observation. Therefore, the dominant contribution to C probably has to be generated after infiation. 

Finally, we mention that in Ref. [7§|, more general vector infiation models are constructed, with the mass term 
replaced by a more general potential. These models are again equivalent to models with a large number of scalar 



With a standard cosmology after inflation, this high inflation scale indeed corresponds to N ~ 55. 
#26 'pijis means that /^matter is obtained from a fiat spacetime expression using the equivalence principle even though the dependence of 
/ upon the bosonic fleds means that the full lagrangian violates the equivalence principle. Using the vierbein formalism, £mattor can 
include spinor fields. 
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fields. The spectra of the field perturbations are the same as before (since they invoke only almost exponential infiation 
without specifying its origin) but their effect on depends in general on what happens at the end of inflation [45| . 
which is determined by other terms in the action. 

X. CONCLUSIONS 

Until recently, it has been assumed that only scalar flelds play a signiflcant role during inflation. Then the spectrum 
of the curvature perturbation is statistically isotropic and homogeneous, and so are higher correlators that would 
correspond to non-gaussianity. Now, it is being recognised that vector fields might be significant during infiation. In 
that case, the correlators of the curvature perturbation will at some level be anisotropic (though still homogeneous). 
The anisotropy will occur if an unperturbed vector field causes anisotropy in the expansion rate, because that will 
cause the correlators of the scalar field perturbations to be anisotropic. It will also occur if a vector field perturbation 
contributes significantly to the curvature perturbation. 

In this paper we have for the first time given expressions for the spectrum and bispectrum of the curvature pertur- 
bation, which include the second of these effects for a generic vector field. 

On the theoretical side, we have for the first time considered the generation from the vacuum of a longitudinal 
vector field component, which will be present in the absence of gauge invariance. Taking its action to be that in 
Eq. (|78l) , we have shown that it can be described by a quantum field theory, according to which its spectrum is twice 
that of the transverse field components. 

We have also given general formulas for the statistical anisotropy of the spectrum and bispectrum, in terms of the 
longitudinal and transverse spectra of the nearly-gaussian vector fields. On the observational side, this leads to a very 
interesting situation regarding statistical anisotropy, which is very similar to that obtained a few years ago regarding 
non-gaussianity. The accepted mechanism for generating C, from the perturbation of the field(s) responsible for slow 
roll infiation, predicted negligible non-gaussianity [i^l, and gaussianity was taken for granted in most early analysis 
of the observations. Starting with the curvaton model ji^, [43, it was found j49| that instead the non-gaussianity 
could be large, and this motivated an intensive search for non-gaussianity. 

Now that vector field contributions to the curvature perturbation are under consideration, statistical isotropy, 
which previously was taken for granted, should be reconsidered. We look forward to the opening up of a new area of 
research, in which predictions for the anisotropy are developed, and confronted with observation. In this context it 
should be emphasised that the bispectrum (and higher correlators) of the curvature perturbation might be completely 
anisotropic*^^, corresponding to the dominance by one or a few vector fields. 
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